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Abstract 

In the Minimal Super symmetric Standard Model with complex parameters 
(cMSSM) we calculate higher order corrections to the Higgs boson sector in the 
Feynman-diagrammatic approach using the on-shell renormalization scheme. The 
application of this approach to the cMSSM, being complementary to existing ap- 
proaches, is analyzed in detail. Numerical examples for the leading fermionic cor- 
rections, including the leading two-loop effects, are presented. Numerical agreement 
within 10% with other approaches is found for small and moderate mixing in the 
scalar top sector. The leading fermionic corrections, supplemented by the full log- 
arithmic one-loop and the leading two-loop contributions are implemented into the 
public Fortran code FeynHiggsFastC. 
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1 Introduction 



The search for the hghtest Higgs boson is a crucial test of Supersymmetry (SUSY) which 
can be performed with the present and the next generation of accelerators. The prediction 
of a relatively light Higgs boson is common to all supersymmetric models whose couplings 
remain in the perturbative regime up to a very high energy scale A precise prediction 
for the mass of the lightest Higgs boson and its couplings to other particles in terms of the 
relevant SUSY parameters is necessary in order to determine the discovery and exclusion 
potential of LEP2 and the upgraded Tevatron, and for physics at the LHC and future 
linear colliders, where eventually a high-precision measurement of the properties of the 
Higgs boson might be possible 0. 

The case of the Higgs sector in the CP-conserving MSSM has been tackled up to 
the two-loop level by different methods such as the Effective Potential (EP) method 
the renormalization group (RG) improved one-loop EP approach and the Feynman- 
diagrammatic (FD) method using the on-shell renormalization scheme The applica- 

tion of different methods lead to thorough comparisons between the different approaches. 
Most prominently the comparison between the RG improved one-loop EP result and the 
FD result and most recently between the FD and the EP result IP,!!^], have been 

performed. These comparisons, showing agreement where expected, lead to deeper insight 
into the radiative corrections in the MSSM Higgs sector and thus to the confidence that 
the higher-order contribution, although being large, are under control. 

In the case of the MSSM with complex parameters (cMSSM) the higher order correc- 



tions have yet been restricted, after the first more general investigations to evaluations 



in the EP approach [^,|TB[ and to the RG improved one- loop EP method |n, r3|. While 
in the MSSM without complex parameters the FD calculation, using the on-shell renor- 
malization scheme, has provided the only complete calculation at the one-loop level 
and furthermore the relevant logarithmic and non-logarithmic corrections at the two-loop 
level a corresponding calculation in the cMSSM has been missing so far. 

This paper provides the next step into this direction: it is shown in detail how the 
FD method, employing the on-shell renormalization scheme, can be applied to the Higgs 
sector of the cMSSM. The general analysis is exemplified at the leading fermionic one-loop 
corrections, showing the applicability of the method and providing the full corresponding 
analytical result. For numerical examples and the comparison with existing approaches, 
the result is supplemented by non-leading corrections at the one- and two-loop level taken 
over from the real MSSM case. All results are finally incorporated into a public Fortran 
code. A detailed analysis, including a full one-loop calculation and the dominant two-loop 



corrections to the cMSSM Higgs sector will be presented elsewhere [|T7 



The rest of the paper is organized as follows. In Section 2 we review the Higgs sector 
and the scalar quark sector of the cMSSM, providing all relevant information about the 
relations of physical and unphysical parameters, the masses and the mixing angles. The 
renormalization in the on-shell scheme in the cMSSM Higgs sector is presented in detail in 
Section 3, together with the analytical result for the leading fermionic corrections obtained 
in this approach. Section 4 briefly reviews the evaluation of the Higgs boson masses and 
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couplings in the FD approach. Numerical results for the comparison with other approaches 
are given in Section 5. Section 6 contains the description of the corresponding Fortran 
code FeynHiggsFastC. The conclusions can be found in Section 7. 



2 Calculational basis 

2.1 The tree-level Higgs sector of the cMSSM 



The (c)MSSM Higgs potential reads |T8 



+ ^^^^ (^1^1 - n^n^f + ^ \n,n2\\ (i) 

o Z 

where m1,ml,ml2 are soft SUSY-breaking terms, g,g' are the SU{2) and f/(l) gauge 
couplings, and ei2 = —1. The doublet fields Tii and I-L2 are decomposed in the following 
way: 



Hi 



Hl\_( v, + {<P\ + ix\)lV2\ 



^ is a possible new phase between the two Higgs doublets. From the unphysical parameters 
in eq. (|l|) the transition to the physical parameters (including the tadpoles) is performed 
by the following substitution (see also Refs. |]6|,pl|): 
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tan/3 is the ratio of the two vacuum expectation values, tan/3 = V2/V1, 



sin/?,C/3 = cos (3, cw = M^/Mz, s 



w 



Ml. 



Ml. 



M^, 



and Sf3 = 
where (as 

will be shown below) Mh± is the mass of the charged Higgs boson H^. Contrary to the 
real case, where the mass of the CP-odd Higgs boson, M4, is used as input parameter, 
in the cMSSM Mh± is chosen as physical parameter, since the field A = spXi + C/3X2 (as 
will be shown later) mixes with the fields 0i and 02- ^1 and t2 denote the tadpoles of the 
fields 01 and 02, whereas Ia is the tadpole of the field A. The expressions for the tadpoles 
can be obtained directly by expanding the Higgs potential eq. (|l|) in the fields from the 
terms linear in 0i, 02 and A. 

In the cMSSM all neutral Higgs bosons can mix. Therefore the following (4 x 4) mass 
matrix has to be considered |TTI|: 



M, 



Higgs 



Ms Msp 
Mg^P Mp 



(4) 



resulting in the Lagrange density 
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(01,02,Xl,X2) M] 



Higgs 
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\X2 J 



(5) 



Here Ms denotes the (2 x 2) mass matrix of the fields 0i and 02 (the CP-even mass matrix 
in the real MSSM), Mp represents the (2 x 2) mass matrix of the fields xi and X2 (the 
CP-odd mass matrix in the real MSSM). Msp denotes the mixing terms (which are always 
zero in the real MSSM). The three matrices are given in terms of physical parameters by 
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(7) 
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with tx = tx e/{2swMw), x = 1,2, A. 

Similarly the matrix of the charged Higgs bosons is given by 



-tiCpsj + t2Sp{l + cj 



2.2 Rotation with /3 

The angle (3 diagonalizes (up to tadpole contributions) the matrix Mp: 



G 
A 



Xi 

X2 



Sp 0/3 



X2 



(Xi,X2) Mp 



Xi 

X2 



(Xi,X2) D{P)D+{P) Mp D{P)D+{P) 1^ 



= {G,A)M^ 



a 

A 



with 



j^D ^ I Cpti + Sph Spti- Cfit2 



This also affects the matrix Msp. Defining the (4 x 4) matrix 

' 1 \ 



D\(3) = 



D{P) i ' 



the rotation of Mniggs can be performed: 
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V X2 / 
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with 

-( ^sp\ nei 



and 



The angle /3 diagonahzes (up to tadpole contributions) also the matrix Mq: 



1 



with 



{G-,H-)MS %\ (19) 



C^tl + S/3t2 +^'3tl-C^t2 \ (20) 



2.3 Rotation with a 

The angle a is defined as 

_ I -| j2 

tan2a = tan2/3^^#^ (21) 

It diagonahzes (up to tadpole contributions) the matrix Ms {s^ — sin a ,Ca — cos a ): 



^^^^ 



(0.02) Ms j 

(0i,02)i^(«)i?+(«) Ms D{a)D+{a) j 
{H,h)Mi(^^^^ (23) 
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with 



Mi 



M 



S^Cf^isl - Cl) + SaCe,{c} - S}) 



+Ml 



+t2 



-C/3S^s2 + 2SaCaS^p + C/jC^ (1 + S 
zlc^Sf^ + Sfsslil + c}) 



Sl3Cp{sl - cl) + SaCa{c} - S^) 



(24) 
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-2c„s„4 + cp{-clsl + sl{l + 4)) 



+ SaCaSpil + 2c|) 



Using the eq. ( pi]) and setting the tadpoles to zero one obtains: 



M 



D 







- CaSpY 

2 ^ 



The rotation with a also affects the matrix Mgp. Defining the (4 x 4) matrix 

D\a) = 

the rotation of M^^ggg can be performed: 



D{a) 
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(27) 



(28) 
(29) 



2.4 Tree-level expressions 

At tree-level all tadpoles can be set to zero. In the 0i-02 sector this ensures that f 1^2 are 
the vacuum expectation values. In the Xi"X2 sector this corresponds to a redefinition of 
the phase of so that the phase e*^ is absorbed . 



One arrives at the following masses at tree-level: 
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(30) 



The entries for the Goldstone bosons G and G^ are to be understood in the Feynman 
gauge. At tree-level there is no CV violation in the cMSSM Higgs sector. The fields h 
and H are decoupled from the fields A and G. 



2.5 The scalar quark sector in the cMSSM 

The mass matrix of two squarks of the same fiavor, and qr, is given by 



M, 



Ml + ml 
mq Xq Ml + ml 



(31) 



with 



Ml 
Ml 



M| + M|cos2/3 {Il-Qqslr, 
M?, + M| cos 2/3 Qgsl 
;U*{cot [3 , tan/3}, 



(32) 



A- 



where {cot /3 ,tan/3} applies for {up, down}-type squarks respectively. In an isodoublet 
the SU{2) symmetry enforces that Mq has to be chosen equal for both squark types. The 
Mq, on the other hand can be chosen independently for every squark type. In the scalar 
quark sector of the cMSSM Ng + 1 phases are present, one for each Ag and one for fi, i.e. 
Ng + 1 new parameters appear. As an abbreviation it will be used 



0g = arg (Xg) . 

As an independent parameter one can trade arg (Ag) = (pAq for 
The squark mass eigenstates are obtained by the rotation 



(33) 





(34) 
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with 

= ( :.!]=( :-../.|:;| ) . s-^' = > . (35) 

where the matrix with 0^ ^ diagonahzes Mq . The mass eigenvalues are given by 

Xq^\Xq 

2 2 I l 

^51,2 = "^9 + 2 



M^ + M^^ J(Ml-Ml)2 + 4m2|X,P , (36) 



independent of the phase of Xg. The unrotated squark mass matrix can now be expressed 
in terms of the physical parameters mq^^rriq^ and the q mixing angle: 

^_^( Cqc*qml + Sqslml 4c*-(m|-m|) \ 

3 Calculation of the renormalized self-energies 

3.1 Renormalization 

The renormalization is performed as follows: 

Ml± Mfj±+SMfj± 

M| ^ M| + 5M| 
tx tx ~\~ Stx, X 1, 2, j4 
tan P tan P + 6 tan /3 

^2 ^ (38) 

The counterterm for the A tadpole can be understood as the effect of re normalizing the 
phase ^ of 0.2 ■ 

In the following we will concentrate on the contributions that are relevant for the 
leading mf corrections (or any corrections of the type ~ m^) for the masses of the neutral 
Higgs bosons. There, only 5Mjj± and 5tx,x = 1,2,^4, enter (see also Ref. [§). 
The renormalized self-energy is then given by 

Re ±H± (0) = (0) - 5M^± , (39) 

the renormalized tadpoles are given by 

ix = Tx + 6tx, x = l,2,A, (40) 

Tx represents the one-loop contribution to tx. 
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The on-shell renormalization conditions are imposed: 

ReS//±(0) = 0, (41) 
4 = 0. (42) 

This results in the on-shell renormalization constants 

5Ml^ = Sj,±(0), (43) 
5t, = -T, . (44) 

Since the charged Higgs boson is renormalized on-shell, its mass does not receive higher- 
order corrections. 

3.2 Renormalized self-energies 



With the on-shell renormalization constants derived in Sect. |3.1| the renormalized neutral 
Higgs boson self-energies read: 



+Ti7TT^^i-^c^s^4 + C/3(-44 + + 4))) 
e 



2swMw 

+7^2 ^^J^^ (-2c.3.4 + 5/3(4(1 + 4) - sl4)) (45) 

r2 



e 

Lw 

+T2— — — (2s„c„4 - clclsp + (1 + cl)slsp) (46) 

thH{^) = J:hH{0)-5M]j4spCp{sl-cl) + SaCaicl-sl)) 

+^^ 9 \j (4(4 - 4) - s^c^c^il + 2sl)) 

+T2-^— (4(c^ - si) + s^c^s^il + 24)) (47) 

±aa{0) = Eaa{0)-6MI^ (48) 
Sgg(O) = ^GGiO) + -^-—{-T,Cf,-T,sp) (49) 

Sag(O) = ^ag{0)+ \ {-T,sp + T^cp) (50) 

6 

S,,a(0) = S/ia(0)+Ta7- 7-^(-CaS/3 SaC/j) (51) 

6 

Sha(O) = SHyl(0)+TA7r T^{-SaS(j - CaCp) (52) 

6 

S/,g(0) = S/jg(0)+Ta- 7-^(-SaS/3 - C„C/3) (53) 

Sj^g(O) = Shg(0)+Ta _ (-gaC/3 + Cag;3) (54) 
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+Ti ^ {-cpslsl + 2s„c„4 + 4c/3(l + 4)) 

e 



3.3 Evaluation of mf contributions 

For the evaluation of the leading mf corrections in the Feynman-diagrammatic (FD) 
approach the diagrams shown in Fig. |I]have to be evaluated for all self-energies Sst(O), st = 
hh, HH, hH, AA, GG, AG, hA, HA, hG, HG. Concerning the tadpole contributions the 
diagrams of Fig. |^ have to be considered. 




Fi gure 1: Generic Feynman diagrams for the mf contributions to Higgs self-energies. 




Figure 2: Generic Feynman diagrams for the mf contributions to Higgs tadpoles. 

For sake of simplicity we now switch back to the 0i-02 basis (i.e. a = 0), where the 
results have a much simpler form. The corresponding results in the h-H basis can be 
obtained by the rotation 



T^hh = sin^ a + cos^ a — 2 sina COS a S(^j02 

^HH = cos^ a t,(f,-^ + sin^ a S<^2 + 2 sin a cos a S0i</,2 

TjhH = —sin a cos a (Tj^^ — Tifi,^^ + {cos^ a — sin^ a )Ti^^^2 

^hA = - sin a S^-^A + cos a J^^^a 

t^HA = cos a S^-^A + sin a t^^^A 

^hG = -smat^^G + cosat^^G 

Ehg = cosa S^^G + sina S<^2G'- (55) 
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The corrections have been obtained using the program FeynArts 3 [1^, employing 
the recently completed MSSM model file ||20|^. Details about the calculations with Feyn- 
Arts can be found in Ref. [^. In the approximation of the leading corrections (and 
with rub = 0) the result for the renormalized self-energies of eq. (^) - 



reads: 
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(by renormalization) 
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^ Only the non-SM counterterms had to be added. 
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(58) 

(59) 

(60) 
(61) 
(62) 
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(67) 
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As expected, S^g(O) = S/fG(0) = Sag(O) = Sgg(O) 
decouples [0. In order to show the finiteness of S^t, st ■ 



0. 



i.e. the Goldstone boson G 

J. J. ^ ^ ^2, 01 A 02^ 



it was necessary to employ the SU{2) symmetry in the scalar quark sector, see Sect. |275 . 
In the simplified case of the leading mf corrections and = (i.e. no mixing in the 
b sector) only the left-handed scalar bottom quark, contributes, where its mass is given 
by eq. (^). 

Exactly analogous expressions have been obtained for the leading corrections (with 
mi m^ and S/3 cp (except in the Ah±a prefactor)), which can be relevant for large 
tan/3. Analogous to the mf corrections also the corresponding contributions (up to 

have been evaluated. 



the color factor and with 



m^ 



A main difference compared to the RG improved EP approach as presented in Ref . 1 14 1 
is the validity of the result as a function of the t sector parameters. Since the FD result is 
obtained directly in terms of the physical parameters in the squark sector, the results of 
the FD approach are valid for arbitrary mixing in the t sector, whereas the RG method 
is restricted to (m| — m?J/(m| + m|) ^1/2. 



3.4 Corrections beyond one-loop order 

Since it is known in the case of vanishing complex phases that the two-loop corrections to 
the neutral Higgs boson masses can be large, for the further numerical examples and com- 
parisons as presented in Sect. ^ the leading contributions at OiGpO-s'fni) ^"^^ C'(G|.mf) 
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are taken into account. For sake of simplicity, up to now the two-loop corrections are taken 
over from the CP-conserving case. The leading corrections then only affect S^j^jlO) 
are valid for arbitrary Higgs sector parameters. They are given by P,p2|,p3[ 



''4'24>2 



(0) 



X 



GfV2 




6 



sin P 




777.- ^ 

*^ sin^ 29i 



2 2 

rriz — mj 



+ \ , *^ sin^ 2^,- log 

2ml 




Ms has to be chosen according to 



(69) 
(70) 



Ml, 



Mr 



Mf^Ml + ml {Ml +MfJ+ mj 



(71) 



and rfit denotes the running top quark mass, Wit = mt{mt). Mi^,Mi^ correspond to 
Mq,Mq, in eq. (^) respectively. Contrary to the presented one-loop result for S^^t^j, 
eq. ( |69D is valid only for not too large mass splitting between the two t mass eigenstates, 
but still gives a rather good approximation for a large part of the MSSM parameter 
space ||2^. The full result in Ref. fl^, however, will be obtained in terms of the physical 
parameters and thus be valid for arbitrary mixing in the t sector. Also eq. ([T0|) is valid for 
not too large mass splitting in the t sector p|,|23[| . However, since the numerical effect of the 
correction in eq. ( [70D is at the ~ 2 GeV level 0, this additional uncertainty is neglected. 
Furthermore, eq. ( ffOD has been obtained in the MS scheme, while all other corrections in 
this paper are evaluated in the on-shell scheme. The corresponding uncertainty is only of 
0{a'^as) and expected to be below ~ 1 GeV and therefore neglected. 



4 The neutral MSSM Higgs sector 



4.1 The Higgs boson masses 

In this section, for sake of completeness, we review the derivation of the Higgs boson 
masses from the calculated higher-order Higgs boson self-energies. Since in the approx- 



imation used in Sect. the external momentum has been set to zero, this step of the 
evaluation is equal to the EP approach |[T^-[14|. In the full FD calculation ||1^ the mo- 
mentum dependence, however, is included, which can lead to corrections of 1 — 2 GeV. 
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Since the Goldstone boson G decouples, see Sect. p73| , the fields 0i, 02 and A form 
a closed subspace that can be evaluated on its own. The masses at higher order can be 
obtained from the diagonalization of the matrix 

/ Mu Mi2 Mi3 
Mcr = M21 M22 M23 
V M31 M32 M33 

The diagonalization is performed with the help of the (3 x 3) orthogonal matrix D^: 

( ^ 

(A 01, 02) Mcv 01 
V 02 

= {A, 01, 02) D^D^^^ Mcv D^D^^ 



(72) 





with 
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D 
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(73) 



(74) 



The numerical evaluation of M^-p and has been presented e.g. in Ref. and is 
also listed here for completeness. The eigenvalues of Mcv are given by 



with 



ei 
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— r + 2 
3 



-p/3 cos 



3 

--r + 2^ -p/3 cos (I + Y 

1 / — f(p 2tt 

— -r + 2d —p/3 cos — — 



(75) 



and 



P 



3s -r^ 



r = -Tr {Mcv) 



2^3 
27 



V9 = arccos 



Tr^ {Mcv) - Tr (m^V) 
14 



t = -Det (Mcp) 



(76) 



(77) 



The rotation matrix can be obtained as 



where 



|Xi|/Ai X2/A2 Xg/Ag 

yi/Ai I1/2I/A2 Z/3/A3 
Zi/Ai 2;2/A2 Izsl/Aa 



A,- 



y2 

X2 

1/1 
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Zl 
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Det 
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Det 
Det 
Det 
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M 



33 



M 
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M33 - m 



M 



Mi3 

M33 - m' 



H2 



Ml 
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M 



22 



m 



M 



23 



M33 - ml 



'33 

Mi3 
M23 



H3 



M22 

M12 
M32 

Mi3 

M23 

M21 
M31 

2 



Mil - mjj^ 
M21 



M21 M22-m- 



M 



31 



M 



^^3 



32 



M12 Mii-ml^^ 
M32 M31 



J 

X sign(y2) 
X sign(z3) 
X sign(xi) 
X sign(z3) 
X sign(a;i) 
X sign (1/2) 



4.2 The Higgs boson couplings 



(75 



(79) 



Again we follow the prescriptions as given in Ref. fl^. Taking complex phases into ac- 
count, all three neutral Higgs bosons are composed of a CP-even part, thus all three Higgs 
bosons can couple to two gauge boson, VV = ZZ, W^W~ . The coupling normalized to 
the SM value is given by 

gH^vv = cpDl^_^ + spDl^_^ . (80) 
The coupling of two Higgs bosons to a Z boson, normalized to the SM value, is given by 



The Bose symmetry that forbids any anti-symmetric derivative coupling of a vector par- 
ticle to two identical real scalar fields is respected, gUiHiV = 0. 
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Finally the decay width of the Hi to SM fermions can be obtained from the decay 
width of the SM Higgs boson by multiplying it with 

{akffy + ialffy], (82) 

with 

9%,un = Dl^-Jsp. gl^uu = Dl^-i C(3/Sf3 (83) 

9mdd = Dl^-i/ci3, 9^H,dd = Dii-iSplcp (84) 
for up- and down-type quarks respectively. 



4.3 The special case of vanishing phases 



In the CV conserving case, e.g. for the leading corrections (pt = <Pn = 0, the CP-even 
Higgs bosons (denoted as h and H with rrih < mn) and CP-odd Higgs boson (denoted 
as A) do not mix. The unrotated mass matrix is then given by 



M, 



CV 



( Ml^ 



V 
















ii - S0,<^i(O) 
(0) 



"^12-^02.^2(0) j 



^5) 



where the square of the CP-odd Higgs boson mass is given by M\ = M^± — S^a(O). For 
a large part of the MSSM parameter space the mass ordering for the three Higgs boson 
masses is given as mn > Ma > mt, i.e. 






Ml 







ml 



and 




(86) 



The mass ordering in eq. ( |7^ can thus imply that in the limit of vanishing phases H2 is 
the CP*-odd Higgs boson. 



5 Numerical examples and comparison with other 
approaches 

The results obtained in Sect. |37^ , eqs. (|56|) - (|67D, have been compared analytically with 
the corresponding results presented in Ref. [jl3l (eqs. (11) - (18c)). Ref. [jl3l calculates the 



leading corrections to the Higgs boson mass matrix in the EP approach. In the approx- 
imation of zero external momentum as applied in Sect. p.3| , the leading mf corrections 
as presented in eqs. (|^) - (pT]) should therefore agree with the corresponding results in 
Ref. [0. Differences due to different renormalization schemes are only expected from 
two-loop order on, see Ref. MM- Complete analytical agreement between the two results 
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is found, if the correction to Ain ,the A boson propagator in eq. (11) of Ref. [|T3l is 

aa ^ ^ 

identified with our renormahzed A boson self-energy, I^aa, given in eq. (|55D. T,aa exhibits 
an additional term compared to the correction to Ain , arising from the fact that in 



Ref. [0 the charged Higgs boson sector has been neglected, while in our approach Mh± 
is chosen as an input parameter, thus introducing T,h± into the result. Therefore, this 
difference only reflects the fact of a different choice of input parameters. A similar ob- 
servation has already been made in Ref. |jl3|, while comparing with Ref. [Q (where also 



analytical agreement in the appropriate limits has been found.) 

In the following subsections some numerical examples are presented and compared 
to results obtained in the RG improved EP calculation. The examples are based on 
the results given in Sects. |3]^, |3]^. They are meant to illustrate the possible effects of 
complex phases in the MSSM. For a phenomenological analysis, however constraints on 
CP-violating parameters from experimental bounds on electric dipole moments (EDMs) 
have to be taken into account, see Sect. |^. On the other hand, the bounds from EDMs 
can easily be evaded by making the first two generations sufficiently heavy ||2^. A more 



detailed phenomenological analysis of the FD results, including the full one-loop and 
leading two-loop corrections in the cMSSM to the Higgs boson self-energies, and taking 



into account all existing experimental constraints can be found in Ref. [|T7 



5.1 Higgs boson masses 

In Fig. 1^ the two lightest Higgs boson masses, rriHi ^"^^ ^H2i shown as a function of the 
phase of the trilinear coupling in the t sector, ipAf The soft SUSY-breaking parameters 
are chosen to emphasize the effect of the CP- violating phases, Msusy = 500 GeV, = 
1000 GeV and = 2000 GeV. The phase of fi is chosen to be zero, except for the 
lower right plot, where it is set to ip^ = tt/2. The phases in the h and r sector are set 
to zero. The different plots show the variation with tan/3, tan/3 = 2,5,20. In the CP- 
conserving case for the above chosen soft SUSY-breaking parameters, tan /? = 2 is already 



excluded by Higgs boson searches ||25|. However, in the CP- violating case this limit 
is weakened due to possible suppressed production cross section and/or suppressed 
decays of the lightest Higgs boson to h quarks, see Sect. |5]^. In each plot different values 
for the charged Higgs boson masses have been chosen, Mh± = 150, 200, 300, 500 GeV. 
The largest effects of the phases are observed for small tan f3 and small Mh± . For large 
Mh± the effects of the CP-violating phases become negligible small. 



A numerical comparison with e.g. Fig. 3 in Ref. [|1^ shows agreement better than 



10% for not too large phases, (pAt ^ 0.8. A larger phase corresponds to larger mixing 
in the t sector. This, on the one hand, makes the corrections and thus the uncertainties 
in the Higgs sector larger. On the other hand, the RG improved EP calculations tends 
to loose accuracy for too large mixing in the t sector. The agreement improves slightly 
if a comparison with the more complete result of Ref. |jl5| (see e.g. Fig. 1) is performed. 
Furthermore, it has been shown in Ref. 0] that differences in the Higgs boson masses 
arising from different renormalizations can be significant, especially for large t mixing. 
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Therefore agreement better than 5-10% cannot yet be expected for all parameter sets 
due to the different renormalizations employed and the yet more complete evaluation 
performed in the RG improved EP calculation. 

In Fig. I the mass difference of the two heavier Higgs bosons, niH^ — "^^2 1 is shown as 
a function of (^Af The other parameters are chosen as in Fig. ^. A large enhancement of 
the mass difference can be observed for small tan/3. The agreement with Refs. |14, 15] is 
found at the same level as for Fig. 0. 



5.2 Higgs boson couplings 

In Fig. H the coupling of the lightest Higgs boson to two SM gauge bosons, relative to its 
SM value, is shown as a function ^p^f The other parameters are chosen as in Fig. ^. Large 
suppressions occur for small values of Mh±- For Mh± ^ 250 GeV no suppression could 
be observed. For small tan/5 the suppression can amount several orders of magnitude, 
whereas for large tan (3 a suppression by a factor of 10 can be observed. These results can 
be compared with the RG improved EP approach, Refs. Fig. 5 and Ref. |]T3[ Fig. 1. 
As for the Higgs boson masses, we find reasonable agreement for not too large values of 

In Fig. H the decay rate of the lightest Higgs boson to h quarks, V{Hi bb), relative 
to its SM value, is shown as a function (pAt- The other parameters are chosen as in Fig. |^. 
The MSSM decay rate, although dependent on the complex phases, is considerably larger 
than the SM one for most parts of the parameter space. This renders the bb channel the 
main decay channel also in the cMSSM. 



6 The Fortran code FeynHiggsFastC 

The results presented Sect. and Sect. are incorporated into the Fortran code Feyn- 
HiggsFastC. They are supplemented by the subleading one-loop corrections from the 
t/t sector ||22[ as well as by the full logarithmic one-loop corrections from all other sectors 
of the MSSM, obtained in the RG approximation [|]. 

In the front-end of the code, the user can specify the input parameters, including all 
relevant complex phases. This part can be manipulated at the user's will. The main 
part of the code consists of the routines needed for the evaluation of the higher-order 
corrections to the neutral Higgs boson mass matrix, and should not be manipulated. 
FeynHiggsFastC evaluates the following items in the cMSSM Higgs sector: 

• the three neutral Higgs boson masses 

• the effective couplings of one neutral Higgs boson to two SM gauge bosons and of 
two neutral Higgs bosons to a Z boson 

• the changes in the branching ratio for a Higgs decaying to SM fermions 
Furthermore the following "check items" are evaluated: 
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Figure 3: The two lightest neutral Higgs boson masses are shown as a function of (pAt 
for different values of Mfj± . In the first three plots (pfj, is set to zero and tan /3 is chosen 
as 2, 5, 20. In the last plot ip^ = ti/2 and tan /9 = 20 is taken. The other parameters are 
MsusY = 500 GeV, |A| = 1000 GeV and = 2000 GeV. 
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Figure 4: The mass difference between the two heavy Higgs boson, rriHs —'rriH^-, is shown 
as a function of (pAt for different values of Mh± ■ In the first three plots (p^ is set to zero 
and tan/3 is chosen as 2,5,20. In the last plot = 7r/2 and tan/3 = 20 is taken. The 
other parameters are MgusY = 500 GeV, | A| = 1000 GeV and = 2000 GeV. 
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Figure 5: The coupling of the hghtest Higgs boson to two gauge bosons (relative to its 
SM value) is shown as a function of (pAt for different values of Mh±. In the first three 
plots (ffj^ is set to zero and tan/3 is chosen as 2,5,20. In the last plot 9?^ = 7r/2 and 
tan (3 = 20 is taken. The other parameters are Msusy — 500 GeV, \ At\ — 1000 GeV and 
\n\ = 2000 GeV. 
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Figure 6: The decay rate of the hghtest Higgs boson to b quarks, T{Hi bb), relative 
to its SM value, is shown as a function of (pAt for different values of Mh±- In the first 
three plots is set to zero and tan /3 is chosen as 2, 5, 20. In the last plot tp^ = n/2 and 
tan (3 = 20 is taken. The other parameters are Msusy — 500 GeV, \ At\ — 1000 GeV and 
\n\ = 2000 GeV. 
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the SUSY corrections to the p-parameter, coming from the t/h sector. (The complex 
phases enter only via their effective change of the t and h masses, where they can 
enlarge the splitting and increase the contribution to Ap.) The SUSY corrections are 
implemented in 0{a) and 0{aas)., where the gluino-exchange corrections, which go 



to zero for large rrig have been omitted . A value of Ap outside the experimentally 
preferred region of Ap^^^^ ~ 3 x 10^^ |27] indicates experimentally disfavored t and 



h masses. 



• the EDM of the electron and the neutron, following the calculation of Ref. [p^ p| 
with the convention of common soft SUSY-breaking parameters for up- and down- 
type squarks. Values outside the experimentally allowed ranges indicate either too 
large CP- violating phases or demand heavier squarks in the first two families p4 |. 

The code can be obtained from the FeynHiggs [B^ home page: www.feynhiggs.de . 



7 Conclusions 

We have presented the application of the Feynman-diagrammatic method and the on-shell 
renormalization scheme to radiative corrections in the Higgs sector of the MSSM with 
complex phases. This provides a complementary method to the (renormalization group 
improved) Effective Potential approach that has been used so far for phenomenological 
analyses. The presented set-up can then be used for a detailed study of the cMSSM Higgs 
sector in the FD/on-shell approach. 

The general FD/on-shell method has been analyzed. Details about the renormalization 
in the on-shell scheme and the derivation of the renormalized Higgs boson self-energies 
have been presented. As an example the leading fermionic corrections to the cMSSM 
Higgs sector have been calculated analytically, making use of the recently completed 
MSSM model file for FeynArts 3. After showing the generic applicability of the approach, 
some numerical examples have been calculated. The leading fermionic corrections have 
been supplemented by the leading two-loop corrections. Results have been obtained for 
the masses of the neutral cMSSM Higgs bosons, their couplings to SM gauge bosons and 
their couplings to SM fermions. Reasonable agreement better than 10% with the RG 
improved EP method has been found for not too large mixing in the scalar top sector. 

Finally the public Fortran code FeynHiggsFastC has been presented. It provides the 
evaluation of the masses and couplings of the cMSSM Higgs bosons in dependence of the 
relevant cMSSM parameters, including all possible complex phases. Besides the leading 
fermionic one-loop and the leading two-loop corrections, also the full logarithmic one- 
loop contributions, taken over from the real MSSM, have been implemented. The code is 
obtainable at www.feynhiggs.de . 

^ We thank C. Schappacher for providing the corresponding Fortran code. 
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